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Abstract 

We demonstrate that the distinctive features of the forward differential cross 
section of diffractive leptoproduction of a vector meson can be legitimately cal- 
culated in perturbative QCD in terms of the light-cone qq wave function of the 
vector meson and the gluon distribution of the target. In particular, we calculate 
the Q 2 and nuclear dependence of the diffractive leptoproduction of vector mesons 
and estimate the cross section. The production of longitudinally polarized vector 
mesons by longitudinally polarized virtual photons is predicted to be the dominant 
component, yielding a cross section behaving as Q~ 6 . The nuclear dependence of 
the diffractive cross sections, which follows from a factorization theorem in pertur- 
bative QCD, provides important tests of color transparency as well as constraints 
on the shadowing of the gluon structure functions and the longitudinal structure 
functions of nuclei. 
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1. Introduction 

We shall consider in this paper the small momentum transfer coherent electro- 
production of vector mesons, j*(q) + p — > V(q + A) + {p — A), where the target, p, 
can be either a nucleon or a nucleus, and where the state (p — A) is identical to p 
except for a small momentum transfer. Here V can be any possible vector meson, 
p°, u, 0, J/ip, T. We shall be concerned with the kinematic region where x = Q 2 /s, 
and My/s are small while Q 2 /Aq CD is large. Small x means that large longitudinal 
distances, in the rest frame of the target, are involved. The effective longitudinal 
distance during which the process takes place is large: Az jj-qi [1]; where M 
is the target mass and s = (p + q) 2 . Our interest is in the possibility of applying 
perturbative QCD (PQCD) to the calculation of hard processes characterized by 
large longitudinal distances. We demonstrate an interesting interplay of perturba- 
tive and nonperturbative QCD effects in the region where the coupling constant is 
small but distances are large which leads to a new way to probe light-cone wave 
functions of hadrons. In the case where Az ^> 2R, with R the radius of a nuclear 
target, we predict that interesting color transparency effects will occur in diffrac- 
tive electroproduction of vector mesons. In the last section of the paper, we briefly 
discuss QCD predictions for the nuclear dependence of the diffractive cross sec- 
tions and show how such measurements can provide important constraints on the 
shadowing of the gluon structure functions and of the deep-inelastic longitudinal 
structure functions of nuclei. 

In general, the physics underlying our PQCD calculation can be directly tested 
through the striking nuclear effects predicted in vector meson leptoproduction 
which differ from those that would result within the Glauber approximation. How- 
ever, the thrust of the present paper is not a detailed discussion of nuclear effects, 
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but rather a general description of the <3 2 -dependence of the forward differential 
cross section, an analysis of the rigor of the PQCD calculation, a semi-quantitative 
estimate of the magnitude of the cross section, a prediction as to which polariza- 
tion dominates the large Q 2 reaction, and a calculation of the nuclear dependence 
of electroproduction of vector mesons. 

The main features of our perturbative QCD analysis agree with those obtained 
in the non-perturbative approach proposed by Landshoff and Donnachie [2] and 
discussed in more detail by Cudell [3]. In the DL model, the Pomeron is represented 
by the effective exchange of two non-perturbative gluons coupling via an effective 
constant coupling and the p wavefunction is approximated as a non-relativistic 
vertex where the quark and anti-quark have equal four-momenta. In our analysis, 
the two gluon aspect of the QCD Pomeron emerges automatically at large Q 2 in 
a form directly related to the proton's gluon structure function. In addition we 
treat the relativistic structure of the vector meson generally. We find that the 
leading twist contribution to p leptoproduction is controlled by the p distribution 
amplitude (p(z,Q), the valence qq wavefunction which controls large momentum 
transfer exclusive processes. As in the DL model, PQCD predicts that the domi- 
nant leptoproduction amplitude couples a longitudinal photon to a longitudinally 
polarized vector meson, and the leading cross section <jl(s,Q 2 ) falls as l/Q 6 . The 
agreement of the data with this form is shown in Cudell's paper. [3] 

We also predict that the cross section for production of a transversely polarized 
vector meson will fall as Q~ s . In fact, it is the end-point contributions (which com- 
plicate the analysis of elastic processes) that yield the dominant contribution to the 
cross section for a transversely polarized V. In contrast, we show that end-point 
configurations are unimportant for the diffractive production of longitudinally po- 
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larized vector mesons. 

Our work is also closely related to that of Ryskin [4] who has made detailed 
calculations of J/vp electroproduction in leading-logarithm in PQCD, employing the 
constituent quark model. In this work we focus on the dominant longitudinally 
polarized channels and find that an analogous analysis can be applied to light vector 
mesons. However, if in our formulae we take a nonrelativistic approximation for 
the wave function of the vector meson, then we find that the cross section for the 
production of the J ftp at large Q 2 is a factor of 4 less than that calculated in Ref. [4] 
in the nonrelativistic approximation. Unlike the heavy quarkonium case where the 
amplitude is controlled by the wavefunction at the origin and the lepton pair decay 
constant, light hadron leptoproduction probes the shape of the minimal Fock-state 
wavef unctions, i.e. the hadron distribution amplitudes (f>y(z,Q) as defined in Ref. 
[5]. This dependence allows tests of non-perturbative predictions from QCD sum 
rules and lattice gauge theory. We also find that this sensitivity may help explain 
the pattern of SU (3)-flavor symmetry breaking seen in the leptoproduction data. 

The electroproduction of vector mesons has also been recently discussed by 
Kopeliovich et al. [6] (KNNZ) in the context of their constituent quark model 
approach to high energy, small momentum transfer processes in which the con- 
stituent quarks interact perturbatively with the nucleon. Indeed, there is much in 
common between the KNNZ approach and the present discussion, including their 
use of Eq. (20) for a nucleon target as derived and applied in Refs. [7] and [8]. 
However, the application of the nonrelativistic quark model to hard processes is 
questionable for this application, since in QCD hard processes should be calcu- 
lated through the distribution of current, not constituent quarks [5]. Moreover, we 
find that the application of PQCD is only legitimate for the production of longi- 
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tudinally polarized vector mesons. Thus, in distinction to Ref. [6], we expect a 
different nuclear dependence for the production of transversely and longitudinally 
polarized vector mesons. Further, the eikonal approximation used by KNNZ for 
the nuclear case is at variance with Eqs. (33) and (37), which are derived from 
factorization in PQCD. As a result we expect that nuclear effects will be lead- 
ing twist — logarithmically decreasing at large Q 2 and increasing at small x. In 
contrast, KNNZ suggest, on the basis of their constituent quark model, that the 
nuclear effects in diffractive leptoproduction are a higher twist effect. 

In order to achieve a simple result for the forward differential cross section, see 
e.g. Eq. (34), we find it necessary to work in the leading in 1/x and the leading 
in Q 2 /Aq CD limit. As discussed in the body of the paper we believe that it is also 
possible to calculate da/dt at small t in the leading in 1/x approximation in terms 
of the unintegrated gluon distribution discussed sometime ago by Catani, Ciafaloni 
and Hautmann [9], without using the leading in Q 2 /Aq CD approximation. How- 
ever, such a calculation is not likely to lead to such a simple result as Eq. (34), 
which follows in the leading-double-logarithmic approximation. Thus, we have not 
pursued the single-logarithmic calculation, although it is clearly interesting to do 
so. At the same time we will explain at the end of Section 2.2 that corrections to 
the expression we obtain are numerically small. 

Our final results should be taken with some caution. Since they are based 
on a leading-logarithmic, even leading-double-logarithmic, calculation, the nor- 
malizations may not be completely reliable. However, the Q~ 6 dependence, the 
dominance of longitudinal polarization for both the virtual photon and the pro- 
duced vector meson, and the proportionality of the cross section to the unintegrated 
gluon distribution in the target, and therefore the nontrivial dependence of cross 
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section on atomic number, are firm predictions which should not depend on our 
logarithmic approximations. 



2. The Diffractive Cross Section in QCD Near t=0. 

In this section the near-forward differential cross section for 7* (9) + p — > V(q + 
A) + (p — A) will be calculated in QCD. The target, labelled by its momentum p, is 
scattered into a state of momentum p — A which we assume to be a particle of the 
same species as p. For example, p and p— A may both refer to protons. V is a vector 
meson of mass my. We assume s/rriy ^> 1, s/Q 2 ^> 1 and — £ = — A 2 Q 2 , and 
s — (p + <?) 2 ? <? 2 — — Q 2 - We also suppose Q 2 /Aq CD , and Q 2 /m/y are both much 
greater than one. 

The differential cross section for the process described above is 



where A is the polarization of the virtual photon and A' is the polarization of the 
final state vector meson. In the large s and large Q 2 limit, but with s/Q 2 3> 1, 
we expect the amplitude M. to be dominated by two-gluon exchange, a particular 
graph of which is illustrated in Fig. 1. We will prove this statement for the 
production of longitudinally polarized vector meson. In Fig. 1, lines k' and q — k' 
are the quark and antiquark making up the vector meson. 
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2.1. Setting up the Calculation of the Matrix Element 

In proceeding to calculate M. it is useful to view the process in a physical way. 
To that end, we choose a frame where p is essentially at rest (i.e., p+ <C q+) and 
where 

q=(q + ,q_,q ± ) = ^ + ,Z^!,o) (2.2) 
with q + = go + <?3, q~ = qo — q%. Then (q + A) 2 = m v and (p — A) 2 = M 2 give 

(Q 2 + m 2 +A 2 ) A 2 
A- « ^ ^ ^, A + « (2.3) 

The polarization vectors are 

e>(±) = (0,0,e[), ^) = (^> ± ) (2.4) 
for the virtual photon and 



where we have dropped terms proportional to Aj_ in (2.5). 

The process illustrated in Fig. 1 takes place, sequentially in time, as follows. 

(i) The virtual photon breaks up into a quark-antiquark pair with a lifetime t% 
given by 

9+T " =Q2+ fr=f K Q2 ■ (2 - 6) 

Here m is the current quark mass. This estimate is valid for the production 
of a longitudinally polarized vector meson only. In the case of a transversely 



polarized vector meson, the end-point non-perturbative contribution arises 
from the kinematical region where z is close to or 1: z, 1 — z ~ m 2 /Q 2 as 
in the aligned-jet model of Ref. [10] adapted to QCD in Ref. [11]. (See also 
the discussion below.) 

(ii) The quark- ant iquark pair then scatters off the target proton. 

(Hi) The quark- ant iquark pair then lives a time tj determined by 

, k 2 + m 2 , 

before the final state vector meson is formed. We note that ry > T{. 

Thus, the amplitude M can be written as a product of three factors: (i) the 
wavefunction giving the amplitude for the virtual photon to break into a quark- 
antiquark pair; (ii) the scattering amplitude of the quark-antiquark pair on the 
target; and (iii) the wavefunction giving the amplitude for the scattered quark- 
antiquark pair of flavor / to become a vector meson. Following the conventions of 
Ref. [5], we have 

x u J V**) { { (2.8) 

where Ai and A2 are the helicities of the quark-antiquark pair which are conserved 
during the scattering off the target. ip v and ^ 7 are the light-cone wavefunctions 
in the notation of Ref. [5]. We have explicitly extracted the sum over the iV c 
colors of the quarks and the 1/y/Nl from the color singlet normalization of tp v . 
Thus, ip v and T correspond to the wave function and scattering amplitude for a 
quark-antiquark pair of definite color. 



We only obtain a simple result, in terms of the gluon distribution of the proton, 
in the leading-logarithmic approximation (in longitudinal and transverse momen- 
tum) for T. In this leading-logarithmic approximation, the time of scattering of 
the quark-antiquark pair with the target is much less than t\ so that T is effectively 
given as the on-shell scattering of a quark-antiquark pair off the target. Our task 
now is to evaluate ?/> 7 and T and thus to express the amplitude, -Mf, in terms of 
integrals over the exclusive wavefunction ip v . 

2.2. ^ 7 AND THE Scattering of the Quark- Antiquark Pair by the 
Target 

The evaluation of ?/> 7 at lowest order in e is straightforward. In the convention 
of Ref. [5], 

« A (*±, *) = ee f ^ h f V ^ q - k) , (2-9) 

which leads to 

In Eq. (2.10), e is the charge of the proton and ey is the charge of a quark of flavor 
/, as a fraction of the proton's charge. We omit the label / on ^ 7 for simplicity 
of notation. We see from (2.4) that longitudinal polarization apparently gives the 
dominant contribution, though we must wait until we have calculated ip v to see 
that this is indeed the case. Anticipating this result we set e 7 = e 7 (£) and, using 
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(2.4), obtain 

Q 2 + 



K^>*)= (2 - n) 



z(l-z) 

Now consider the scattering of the quark-antiquark pair by the target. The 
relevant light-cone perturbation theory graphs are shown in Fig. 2. (It is important 
to note our convention of using the momenta q — k! and k' + A that enter ip v to 
define the momenta of the lines.) As we have mentioned before, a simple result 
emerges only in the leading-logarithmic approximation in £n Q 2 /s. In that case, 
£+/k + <C 1 and the dominant couplings of the lines £ and £ + A to the quark- 
antiquark pair occur with a 7+ (in a covariant gauge). Thus, the vertices of the 
lines £ + A and £ with the quark-antiquark pair are exactly the same for each of 
the graphs in Fig. 2. Further, all energy denominators are dominated by the £ and 
£ + A lines so that the energy denominators are also the same for each of the graphs 
in Fig. 2. The differences between the different graphs in Fig. 2 are only in the 
labelling of the momenta on the left-hand side of the diagram as they emerge from 
the ?/> 7 wavefunction. The result obtained in the Aj_ <C £j_ limit being considered 
is very simple: 

Tx^ik't, z ] k±,z) = 

x 16tt 3 J {25(A4 - k±) - 5(k' ± -k± + £±) - S(k' ± -k±- £±) } ^ 12 ) 

where I(£) represents the gluon propagators and zero-angle gluon-nucleon scatter- 
ing amplitude shown in Fig. 1. 

Before attempting to evaluate the d 2 £±d£+ integral involving X, let us first do 
the integrals over d 2 k' ± dz' /IGir 3 indicated in (2.8). Using (2.12) one arrives at the 
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combination 



2^ XiX2 (k' ± , z) - r XlX2 (k' ± + t ± , z) - ^ XiX2 (k' ± - t ± , z) = A^ iA2 (2.13) 
which, using (2.11), gives 



_ -2ee / ^ Al ,-A 2 Q^. 



Ai A2 



when £ 2 JQ 2 < I. Thus (relabelling k' ± -> fcjj, 



2 



(2.14) 



1 

-M/ = V^c £ Al , A2 / d^ft^x,*) X(£) A^ iAa . (2.15) 



If instead of a final state vector meson we were considering a virtual photon 
identical to the initial state photon, the exact same dependence on £±_ and £+ would 
appear. This allows one to identify an integral over I(£) with the gluon distribution. 
Indeed, in the leading- logarithmic approximation in £n 1/x and £n Q 2 /A 2 , where 
X is purely imaginary, 



/ 



d2 ^ M+ l l m = i^IS^ ( 4 + 0»)*G(*,g») , (2.16) 



16ir J ^ w N, 



C 



where the factor of Tr/N c arises by virtue of averaging over the color and matching 
anti-color of the initial quark and antiquark. (In the usual convention Tr = \.) 
The simplest way to determine the normalizing factor in Eq. (2.16) is to compare 
the integral on the left hand side of (2.16) with the known relation of the longitu- 
dinal structure function with the gluon density. Once one is sure that the gluon 
distribution should emerge from the integration on the left-hand side of (2.16), the 
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normalization is most easily set by taking the target to be a quark and calculating 
the lowest order contribution to X. However, to establish that x = Q 2 /s and Q 2 
are the appropriate arguments in xG(x, Q 2 ) requires employing the leading in 1/x 
and in Q 2 /Aq CD approximations. 

We can understand these values from the simple picture on which Eqs. (2.8) 
and (2.15) are based. Namely, we must have a strict sequence of events in which (i) 
the virtual photon breaks up into an approximately on-shell quark- ant iquark pair; 
followed by (ii) the scattering of the on-shell quark-antiquark pair by the target; 
and, finally, (iii) the scattered quark-antiquark pair turns into the vector meson. 
This picture requires that the time of scattering of the quark-antiquark pair must 
be much less than tj and ry defined earlier: i.e. r <C t\ < Tf. Using (recall that 
we employ OFPT here) 

this condition reduces to 

U « ^ • (2.18) 

Meanwhile, the dominance of the imaginary part implies that the — components of 
the four-momenta are approximately conserved from the initial to the intermediate 
state where the p + l line is cut, leading to (p+£)- ~ p_ (neglecting — components 
of order l/q+). With this, we compute s — (p + £) 2 ~ M 2 + p~£ + i^s/Q 2 , 
with the last inequality coming from employing (2.18). In the leading in 1/x 
approximation we allow integration over the parton-gluon scattering subprocess 
energy, s, up to this upper limit. Recalling that the x g argument of G(x g , Q 2 ) is 
set by x g s ~ t\, we see that in the leading In 1/x approximation x g should be 
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identified with x = Q 2 /s. To determine the appropriate argument Q 2 , we note 
that the simple form for A?/> 7 given in (2.14) depends on 

J « 1 ■ (2.19) 

In leading in Q 2 /Aq CD we integrate i 2 ^ up to the maximum allowed by this rela- 
tion, and it is this maximum which determines the argument Q 2 appearing in G. 
The result is obviously Q 2 ~ Q 2 . (We note here that so long as A^/Aq^ <C 1, 
the left-hand side of Eq. (2.16) does not depend on A. From (2.3), we see that 
A_ is comparable to = i\/ i+ only when £ 2 L /£+ reaches its smallest value, 
Q 2 /q+ ~ A_.) It is useful to note that these conclusions match closely those ob- 
tained in Ref. [7] where the scattering of a quark- ant iquark pair off a target is 
given by 

2tt 2 

°(b 2 ) = — {b 2 a s (Q 2 )xG N (x, Q 2 )) - =1/s62 Q2=1/fe2 , (2.20) 

where, for simplicity, the result has been stated for the case x q ~ Xq ~ 1/2. 

Can one do better than the leading double log approximation? The answer 
should be yes. One should be able to eliminate the restriction (2.19) and derive an 
expression for M. in terms of integrals over the unintegrated gluon distribution and 
the vector meson exclusive wavefunction. The result, however, will be significantly 
more complicated than the answer we are about to give for the leading double- 
logarithmic approximation. At the same time, the corrections to (2.16) resulting 
from the elimination of the restriction (2.19) are numerically small. A simple way 
to justify this statement is to perform the calculation in impact parameter space 
where non-£n Q 2 /Aq CD corrections arise from the decomposition of the matrix 
element of the real part of e tej -' bj - — 1 between the wave functions of the virtual 
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photon 7* and a vector meson in impact parameter space. Here bj_ is the inter- 
quark distance in the meson wavefunction and £j_ is the transverse momentum of 
the gluon shown in Fig. 2. The second order term in (£± • b±) leads to Eq. (2.15). 
The corrections in question arise starting at the fourth order in the expansion 
contain a factor of 1/4! and are therefore small. 

2.3. The Final State Vector Meson Wavefunction 
Using (2.16) and (2.14) in (2.15), and (s + Q 2 ) ~ s, one obtains 

— 8isit 2 eefa s Tr f d 2 kj_ f dz _ lV , 

o 

where 

,V a \ at v-\(q~k) v u x (k) v 1 

^A-A( fc ±> ^) = N v — 7= 7 • e — ^ (fc±, 2) • (2.22) 

Vl — 2 V- 2 m ^ 

Here, we have introduced a normalization factor Ny for the particular spin state 

of V being considered. (In the conventions of Ref. [5], Ny = l/\/2 for ev(L).) 

We are now in a position to verify that longitudinal polarizations give the 
dominant contribution. From (2.10) and (2.22) it is straightforward to see that 
longitudinal polarizations give the leading contribution so long as Q is much larger 
than all quark and vector meson masses. In this case, using (2.5) in (2.22) gives 

iPL x (k ± ,z)= -N v ^ v (k^z) (2.23) 

which yields (after performing and setting Tr = 1/2) 

o 
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where 

<P V {Q,z) = J^j,V(kj_,z) (2.25) 

is the distribution amplitude for longitudinally polarized vector mesons with the 
restriction k\ < Q 2 understood for the integration in (2.25). In getting the correct 
normalization for the cross section, it will be important to note that this result 
does not yet include the flavor normalization for the V wavefunction. (The spin 
normalization factor is contained in iVy.) 

In order to obtain an absolute normalization for M.f, it is necessary to relate 
4> v to the experimental observable fy defined by 

(0|^.|V) = ^^, (2.26) 
mye^ 

in terms of which the decay width for V — > e + e~ is given by 



r _ 87ra 2 f v 
3my 

We demonstrate in the Appendix that 



*' V( « )= ^,W (2 - 28) 



(2.27) 



Defining rjy as 



W = 2 fdzWz) (2 ' 29) 



we obtain 



8 isn 2 fyea s r] V 

Mf= ( ' Q) ' ( ' 

The parameter rjy is the effective inverse moment of the vector meson distribu- 
tion amplitude that controls the leading twist contribution to the leptoproduction 
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amplitude. Higher particle number Fock state amplitudes such as the qqg in a 
physical gauge have a suppressed coupling to the small-size quark pair; in order to 
compensate their higher mass dimensions, they must be accompanied by further 
powers of 1/Q. Note that both ey and Ny have cancelled out in relating M.f to 
the experimental observable fy. 

In order to estimate the cross section it is convenient to consider two extreme 
examples for the shape of the p distribution amplitude: 



<j>X = r 3 = —z(l - z) (2.31) 



<j)X = -= 15 — ^(l-^)(l-2^) 2 . (2.32) 
V2N V ^/W C e f 

Here 0i is the asymptotic form of the distribution amplitude, while 02 is of the 
form suggested by Chernyak and Zhitnitsky [12] for pions. The normalizations of 
4>Y and 4>2 are chosen for consistency with Eq. (2.28). A simple calculation gives 
r\y — 3 for 0i and rjy — 5 for 02. We emphasize that in QCD amplitudes for hard 
processes are expressed in terms of the minimal Fock component of the light-cone 
wave function of the meson (see e.g. Eq. (2.21)), not in terms of constituent quark 
model components. This is important when the physics of color screening, which 
is relevant to the transition from the non-perturbative to the perturbative regime, 
is accounted for (c/. the discussion in Section 2.6 regarding the difference between 
the production of transverse and longitudinally polarized vector mesons, and the 
discussion of color transparency in Section 2.8). 
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2.4. The Flavor Dependence of the Cross Section 



Thus far, we have not discussed the flavor dependence of the wavefunction of 
the final state vector meson. If V refers to a neutral p meson, then the relevant 
wavefunction is ^i\ u ^ > —\dd>). This would imply a replacement of e/ in our 
calculations of both Aif and fy by ej — > -^( e u — ^d) — l/v 7 ^- However, since e/ 
cancels out when A4f is expressed in terms of fy, as in Eq. (2.30), the result is 
that we may use Eq. (2.30) for the total amplitude without change provided we 
employ the appropriate value for f p . With our normalization conventions, f p ~ 107 
MeV for the p°. If V refers to a J/ip one can simply use ej = e c = 2/3. Again, 
the explicit value of ej disappears if M. f is expressed in terms of /ju. However, 
for the J/^ the value of rjy is expected to differ substantially from that for the 
p since neither of the wavefunctions given in (2.31) and (2.32) is appropriate; a 
wavefunction having z ~ 1/2 would be more suitable. 

2.5. The Differential Cross Section 

We shall now write the differential cross section for p° production. Modifica- 
tions for other lepto-produced neutral vector mesons are straightforward and in- 
volve choosing a distribution amplitude to replace (2.31) or (2.32). Using Eq. (2.30) 
in Eq. (2.1) gives 



da 
~dt 



( 7 *iV -> VN) = 8 ^fv a EM^(Q)Vy^G(x,Q)] 2 ^ 



It is important to keep in mind that this result gives the differential cross section 
for a longitudinally polarized photon to produce a longitudinally polarized p°; i.e. 
it is not spin-averaged over initial photon states. As an aside, we note that the 
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t-dependence of this diffractive cross section is controlled by the quasi-local two- 
gluon matrix element of the nucleon; in principle, it could have a different fall-off 
than the elastic form factors since the momentum transfer is shared by the two 
gluons. 

It is convenient to rewrite Eq. (2.33) in terms of the leptonic width, Ty using 
Eq. (2.27), 



da 
~dt 



{ffN ^ VN) = ^Ty m y^Q)f[fM)? 



since the coherent sum over the contributing flavors is identical for the diffractive 
amplitude and the decay amplitude. Our prediction for the J/ip leptoproduction 
cross section using Eq. (2.34) is smaller than the one obtained in Ref. [4] by a 
factor of 4, if we assume, as in Ref. [4], a nonrelativistic form ~ 8{z — 1/2) for the 
J/ip meson wave function so that 77^ ~ 2. 

We can also use dispersion relations to determine the real part of the ampli- 
tude. For small enough x and large Q 2 the contribution of the real part is not 
negligible since the effective QCD pomeron intercept is above 1. Including the real 
contribution as a perturbation we get 



^TymyaliQWy | (1 + ij T d-) x G T (x, Q) | 2 
a EM Q 6 N 2 

Since the gluon density rises rapidly at small x in the large Q 2 domain, we 
predict a very substantial rise of the diffractive cross section with energy at large 
Q 2 . For example, for Q 2 ~ 10 GeV 2 , we predict a rise of the diffractive cross section 
at small x by as much as a factor of 100 at HERA energies as compared to the 
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d<J TN _>y N 



dt 



cross section measured at CERN. Obviously this effect would also substantially 
modify the Q 2 — dependence of the cross section at large energies. 

Let us now compare Eq. (2.33) with experimental data for p-meson production 
at large Q 2 and small x. Such data are available from EMC [13] and NMC [14,15] 
and also from the E-665 FNAL experiment [16]. We will use the latest data of 
NMC [15] which were obtained at the highest Q 2 with special attention to removing 
backgrounds due to inelastic processes. We also note that the data of E-665 [16], 
which extend to somewhat smaller Q 2 , are generally consistent with NMC data. 
All of these data confirm the important role of the longitudinal contribution at 
large Q 2 from the measurement of the polarization of the p. If s-channel helicity 
conservation is assumed in the transition 7* — > p, this also determines the ratio of 
ox/or- 

To convert the leptoproduction cross section d 2 a/dQ 2 du from pN to the vir- 
tual photoproduction cross section, we use the standard relation cr^(Q 2 ,u) — 
Ta rN (Q 2 ,u), where T = a EM {y - Q 2 /(2M))/(2ttQ 2 £ 2 (1 - e)) and a rN = 
ot + e&L- Using a tot (^i*N — > pN) as determined by NMC from pD data, the 
slope of the t-dependence of the cross section, b = 4.3 ± 0.6 ± 0.7 GeV~ 2 as 
measured by NMC, and the NMC estimate for o"j nco /j/o"t t = 0.55 ± 0.08 and of 
<j ^*N->o n/' j ?y*n^o n ~ 2-0 at Q 2 = 6 GeV 2 and e = 0.8, we can estimate 



da 



L experiment 
'~t*N^p°N 



14 - 27 nb 
GeV 2 



d,t 



for Q 2 = 10 GeV 2 . (In deriving this range we also assumed that o-lI&t is either the 

o o O o 

same at Q £ = 10 GeV z as at Q z = 6 GcV z or increases linearly above Q = 6 GeV 
as ~ Q 2 .) 



20 



Let us now compare this result with the leading-logarithmic prediction of 
Eq. (2.33). The HMRS DO' parameterization of the gluon distribution in the pro- 
ton [17] gives the value a s (Q)[xG(x,Q)\ = 0.67 at the NMC kinematics x ~ 0.06, 
(5 —10 GeV . The sensitivity of the parameterization to uncalculated higher or- 
der terms and the uncertainty in the evolution scale is illustrated by noting that 
cng^G] = 0.76 if we use \Q 2 instead of Q 2 for the arguments of a s and G. Thus 
we obtain a range of predictions 

dcr^KT ,„o^ (13 - 17) - (36 -47)nb 



dt 



t=o GeV z 



The lower range corresponds to r\ p = 3 assuming the asymptotic form of the pi 
distribution amplitude ~ z(l — z). The upper range corresponds to r\ p = 5 assuming 
that the pi distribution amplitude is similar to the CZ wave function of a pion. 
The distribution amplitude suggested by Chernyak and Zhitnitski [12] from QCD 
sum rules for the pi actually corresponds to a narrower quark distribution than for 
the pion, but it is still broader than the asymptotic form. Thus rfjj^ ~ 3.3 — 3.5. 
If we trust our leading-order estimates, then this value of i] pL leads to 



dt 



16 - 23 nb 
t=o, v % GeV 2 



which is close to the empirical diffractive p leptoproduction cross section. More 
generally, p electroproduction will allow us to test determinations of the meson 
distribution amplitudes from lattice gauge theory or other non-perturbative QCD 
computations. 

Equation (2.34) also allows us to predict the ratio of the yields of various vector 
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mesons at Q z S> m 2 ? : 



«"^ +r - v ^ _ Tvi m Vl rfa 
=o T v 2 m v 2 Vv 2 



R(Vi, v 2 ) - d(7 ( 7£+ jLi/ 2+ T) 



dt 



(2.35) 



It can be seen from Eq. (2.35) that we predict R((p, p) ~ 1.0 x e 2 ((f))/e 2 (p) for 
Vcp/Vp ~ 0.9 suggested by CZ [12]. We also find R(J/4>,p) ~ 1.2 x e J(J/*)/eJ(p) 
for r^j/^, = 2 corresponding to 0jm ~ 5(z — 1/2). Here, we have extracted from 
these two R results the ratios of the effective charge squared, where e 2 (V) = 
1/2,1/18,1/9,4/9 for the p,u,(j),^f, respectively (for standard SU(3) wave func- 
tions), which might naively be expected to determine the i?'s. 

Another interesting feature of the QCD prediction for the production of a 
vector meson at large Q 2 is the universality of the t-dependence of the process — it 
is determined by a universal two-gluon form factor, independent of the vector meson 
type. The upper part of the amplitude corresponding to the transition 7* — > V 
is effectively dipole-like at large Q 2 ; i.e. it should depend weakly on t so long as 
—t Q 2 , implying that the t-dependence of the leptoproduction cross section 
primarily reflects the t-dependence of the gluon - nucleon scattering amplitude. 
The slope of this t-dependence should increase slowly with the incident energy due 
to the Gribov diffusion - shrinkage of the diffractive cone. The data on exclusive 
production of vector mesons at high energies supports this prediction — the large 
Q 2 p- meson leptoproduction cross section has a slope da/dt oc e bt corresponding 
to b = 4.3 GeV -2 [14], which is similar to the small slope b ~ 3.5 — 5.5 GeV -2 
observed for exclusive J/^-meson photoproduction. 
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2.6. Non-Perturbative QCD Effects 

The Q 2 dependence of the diffractive vector meson leptoproduction amplitude 
in the PQCD analysis reflects the overlapping integral between the light cone wave 
functions of 7* and the vector meson. Since the wave function of 7^ is oc z or 
(1 — z) when z or 1 — z vanish, the end-point contributions to the cross section for 
the production of longitudinally polarized vector mesons, arising from Eq. (2.8), 
are small and of order ~ Q 2 /rriy(Q 2 ) 6 . Since the wave function of a vector meson 
bound state should be less singular at z ~ and 1 than the wave function of the 
7*, we conclude that for ox the end-point contribution may be neglected. Thus 
the longitudinal cross section can be safely calculated in terms of PQCD. 

The situation is the opposite in the case of the production of a vector meson 
by a photon with transverse polarization. The wave function of a transversely 
polarized photon is constant at z ~ and 1. As a result, in this case the end- 
point contribution is enhanced. To demonstrate the importance of the end-point 
contribution, let us consider the process where a final state photon is produced 
instead of a vector meson; i. e. the process: 

7* + T -> 7 + T . 
For this reaction, the end-point contribution leads to the cross section 



dt 



1 

Q 



(2.36) 



t=o 

Now note that the wave function of a transversely polarized vector meson is less 
singular than that of the photon as z — > 0, 1. In fact, PQCD predicts that the 
qq component decreases at least as fast as oc z or (1 — z) when z or 1 — z are 
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small. Consequently, the cross section for electroproduction of transversely po- 
larized svector mesons should fall at least as fast as l/Q 8 . (This contribution is 
additionally suppressed by a double logarithmic Sudakov-type form factor.) A sim- 
ilar dependence of the cross section arises from other kinematical regions. Thus 
the experimental investigation of the ratio of longitudinally and transverse po- 
larized vector mesons would help to clarify the relative roles of non-perturbative 
Q CD-end-point contributions and hard physics. 

2.7. DlFFRACTIVE LEPTOPRODUCTION ON NUCLEI 

A key feature of the predictions of PQCD for forward diffractive vector meson 
leptoproduction with Xjlm^x 2Ra is the dominance of small size wavefunction 
configurations. The fact that the integration range of in the vector meson 
light-cone wavefunction in Eq. (2.25) extends to ~ Q 2 implies that the important 
qq configurations coupling to the virtual photon have transverse separations 6j_ ~ 
l/Q. Thus, even in a nuclear target, color screening implies that the coherent qq 
system can only weakly interact, and in leading- logarithmic approximation only 
two gluons in light-cone gauge connect the photon-vector meson system to the 
nucleus, as illustrated in Fig. 2. Thus, as predicted by PQCD color transparency 
[18, 11] the outgoing vector meson in effect suffers no final-state absorption, and 
the nuclear dependence of the 7M — > V A forward amplitude will be identical to 
that for the case where the final state system is a virtual photon, 7* A — > 7* A; i.e. 
it will be close to additive in the nucleon number A. We can also understand this 
remarkable feature of QCD from the space-time arguments given above: the final 
state vector meson is formed from the compact qq pair over a long time tj, and 
thus it does not attain its final physical size and its normal strong interactions until 
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it is well outside the domain of the target nucleus. In fact, much of this physics was 
anticipated before the advent of QCD. The possibility that the outgoing absorption 
of the p in the nucleus would be effectively small in large Q 2 leptoproduction was 
actually first proposed by Yennie in 1975 [19]. The observation that the incoming 
photon has point-like behavior and diminished absorption was discussed in terms 
of a "shrinking photon" by Cheng and Wu [20] and by Bjorken, Kogut, and Soper 
[21]. 

We emphasize that the above reasoning is applicable for a longitudinally po- 
larized vector meson only — the effective transverse size of a produced transversely 
polarized vector meson is considerably larger (although still smaller than for ordi- 
nary hadrons, cf. the discussion in Section 2.6). 

Although the vector meson suffers no final state interactions, the forward am- 
plitude Ai(j*A — > VA) is not strictly additive in nuclear number since the gluon 
distribution itself is shadowed. (This effect is similar to the shadowing of diffrac- 
tive production of high-py jets in the it + A — > 2 jets + A reaction, discussed in 
Ref. [8].) In fact, we see from Eq. (2.33) that 

^( 7 M -> VA) oc a 2 s (Q)[xG A (x,Q)} 2 , (2.37) 

where Ga(x, Q) = A ag ( x >^GN(x, Q) is the gluon distribution in the nucleus. Thus, 
the analysis presented in this paper predicts identical nuclear dependence for the 
forward vector meson diffractive leptoproduction cross sections, the longitudinal 
structure functions F^(x,Q), and the square of the gluon structure functions: 



( 7 M -> VA) 



da 



Ff}(x,Q) 



n 2 



G A( X 'Q) _ A 2a g (x,Q) (2 o R) 

G 2 N (x,Q)- A ■ (2 ' 38) 



(Note that at finite energies one has to interpolate the cross section to the unphysi- 
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cal t — kinematical point.) The nuclear gluon distribution is expected to be more 
strongly shadowed than the nuclear quark structure functions at intermediate Q 2 
because of the larger color charge of the gluon in QCD and thus its stronger inter- 
nuclear interactions. Numerical estimates [22] lead to xGa(x, Qq)/AxGn(x, Qq) ~ 
0.7 - 0.8 (0.4 - 0.5) for A = 12 (200) and x ~ 0.01 - 0.03, a result which 
seems to be supported by the recent FNAL data of E-665 [23] . However, at fixed 
x ~ 0.01 — 0.03, shadowing substantially decreases with Q 2 due to scaling viola- 
tion effects [22], which should lead to an effective increase of transparency for p 
leptoproduction at fixed x with increasing Q. 

A nuclear dependence similar to Eq. (2.38) is also expected for forward diffrac- 
tive T and, possibly, J ftp leptoproduction cross sections even at small Q 2 , although 
in this case the value of x and the evolution scale of G(x, Q) is controlled by My 
rather than Q 2 . Thus both heavy and light diffractive vector meson leptoproduction 
can provide basic information on the nature of the gluon distributions in nuclei. 

The calculations given above are applicable to the near-forward production 
of vector mesons. It should be noted that the physics relevant to the nuclear 
dependence of the leptoproduction cross section will change with increasing t. 
We shall give here a semi-quantitative description of the expected behavior. At 
— tR A /3 <C 1, coherent processes dominate the leptoproduction of vector mesons. 
The nuclear dependence of the diffractive cross section at small t (within the diffrac- 
tive peak) can be estimated by multiplying Eq. (2.38) by the square of electromag- 
netic form factor of nucleus normalized to 1 at t — 0. However, if —tR A /3 3> 1, 
incoherent processes, in which the leading vector meson is accompanied by the 
production of other hadrons from nuclear disintegration, will dominate the cross 
section. The existence of nuclear shadowing implies that gluons at small x cannot 
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be associated with individual nucleons. Thus one can have events where momenta 
— t > 0.1 GeV 2 are transferred to each of several nucleons which subsequently frag- 
ment. The effect can be a slower t-dependence of the cross section and a smaller 
energy transfer per interacting target nucleon than for the scattering off a single 
nucleon. The expected A-dependence is intermediate between that expected for 
shadowing of Ga(x, Q 2 ) and A. 

The recent nuclear target p leptoproduction measurements from the E-665 
experiment [16] appear to indicate onset of the color transparency predicted by 
PQCD for incoherent — > pN(A — 1)' reactions. The onset of this phenomena 
is again for Q 2 ~ a few GeV 2 , the same scale at which Bjorken scaling is observed 
in deep inelastic lepton scattering reactions. Preliminary data [15] from the NMC 
also confirm higher values of the transparency ratio for Q 2 > 3 GeV 2 , observed 
in [16] although the NMC data do not indicate a Q 2 variation of transparency 
in their 10 > Q 2 > 3 GeV 2 range. One needs to be cautious in interpreting 
these data directly in terms of PQCD color transparency of the outgoing p. We 
note that the high-Q 2 NMC data correspond to a range of x where the essential 
longitudinal distances are smaller than the nucleus size. Thus transparency in this 
range of kinematics will reflect to some extent the fact that the virtual photon can 
penetrate deeper into nucleus without interaction. As we have emphasized, the 
nuclear dependence of forward diffractive p leptoproduction which is completely 
coherent on the nucleus can provide a decisive test of color transparency. 
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3. Conclusions 

The analysis of diffractive leptoproduction cross sections presented in this pa- 
per extends the domain of PQCD predictions to a new domain of exclusive hadronic 
reactions. The central focus of this analysis is closely related to the calculations 
of the order a s (Q 2 ) leading twist perturbative contributions to gl{v, Q 2 ) and the 
violation of the Callan-Gross relation. Although the momentum transfer to the 
target is small, the virtuality of the longitudinally polarized photon provides a 
point-like probe of the diffractive 7*iV — > VN process. At high Q 2 the ampli- 
tude factorizes in terms of separable components: the perturbative distribution 
amplitude of the virtual photon, the non-perturbative distribution amplitude of 
the outgoing vector meson system, and a non-perturbative two-gluon matrix ele- 
ment of the target closely related to the gluon structure function. The momentum 
transfer dependence of the diffractive amplitude is thus controlled by a new type 
of non-local two-gluon form factor. We note that since the momentum transfer to 
the nucleon is shared by two gluons, the fall-off can be different from the fall-off of 
elastic electromagnetic form factors. 

The factorization analysis can be extended to the diffractive production of any 
vector meson system of mass M. as long as Q 2 > M 2 . The longitudinal cross 
section always falls as l/Q G at fixed M, t and s. The energy dependence of the 
forward diffractive cross section is also universal, reflecting the behavior of the 
square of the gluon structure function at x ~ Q 2 / s. Related formalisms have also 
been applied to exclusive heavy quarkonium photoproduction [4], exclusive pion 
dissociation to two jets [22], and to two-photon diffractive reactions [24]. 

The existing data for diffractive p leptoproduction appear to be consistent 
even at relatively low Q 2 of a few GeV 2 with the leading logarithmic predictions in 
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both magnitude and in the kinematic dependence. A crucial feature of the PQCD 
predictions, which is clearly evident in the data, is the increasing dominance of 
the longitudinal photon and vector meson production amplitudes with increasing 
Q 2 . More careful studies of the longitudinal to transverse polarization ratios can 
lead to insights into the transition between soft and hard components in QCD 
amplitudes. As we have emphasized, precise measurements of the A dependence 
of the diffractive leptoproduction reactions can lead to new insights into nuclear 
shadowing of the longitudinal structure function. 

In our discussion we have noted only a few of the many empirical tests of PQCD 
possible in diffractive leptoproduction. In principle, the study of these reactions at 
HERA will allow tests of the theory over a huge dynamical range in x and Q 2 . The 
intercomparison of the various vector meson channels can also lead to new tests 
and understanding of the non-perturbative structure of hadronic wavefunctions 
and their flavor-symmetry properties. 
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APPENDIX 



The 7 — > V transition is defined by 

(01 W = • (A.l) 

In the frame defined by Eq. (2.2) (where we shall take Q 2 = —ray), and employing 
the time-ordered formalism outlined in this paper, we find 



V2ef v e~< ■ e v 



ray ( 

AlA; 



(A.2) 



where the amplitude T is trivial, 

T = 16tt 3 5(A4 - fc_L)£(z' - 2) , (A.3) 

and ip 1 and tp v are given by Eqs. (2.10) and (2.22), respectively. In the chosen 
frame, it will be most convenient to isolate fy by considering a longitudinally 
polarized V, with polarization vector as given in Eq. (2.5), while taking tL = 1 
and all other components of e 7 to be 0, implying that 

e y. e v = h v {L) = kt. (A .4) 

Meanwhile, in Eq. (2.10) 7 • e 7 = ^7+ while in Eq. (2.22), for L polarization, 7 • e v 
can be approximated by ^7+e^ = — ^7+my/ ( Z+. Using the fact that 

"Ai(fc) 1 ^A 2 (g~fe) ^A 2 (g~ fc) 1 tt Al (fc) „ r 
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we obtain from Eq. (A. 2) the result (with Q 2 = —m v in (2.10)): 



m 



V 



16tt 3 



my 



With = 2 and using the definition Eq. (2.25), we obtain 



dz 4> V {z) 



fx 



W c e f 2V2N v 
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FIGURE CAPTIONS 

1) A typical two-gluon exchange contribution to the amplitude for 7*(g)p — > 
V(q + A)(p - A). 

2) Light-cone perturbation theory graphs for the scattering of a quark-antiquark 
pair by a colorless target. 
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